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1 Introduction 

In this paper we solve the isoholonomic problem in a homogeneous bundle and apply 
this result to the optimal control problem in holonomic quantum computation. In other 
words, this paper has two purposes; first, we solve a mathematical problem which has 
been unsolved for more than a decade since it was initially proposed by Montgomery ^J. 
Second, we provide a scheme to construct explicitly an optimal controller for arbitrary 
unitary gate in holonomic quantum computation [21 El- 

The isoholonomic problem is one of generalizations of the isoperimetric problem. The 
isoperimetric problem, also known as Dido's problem, is originally proposed in the context 
of plane geometry; what is the shape of a domain with the largest area surrounded by a 
string of a fixed length? The solution is a circle. The isoperimetric problem has a long 
history and various generalizations thereof have been proposed. 

The isoholonomic problem is formulated as follows. Assume that we have a principal 
fiber bundle (P, M, it, G) with a connection. The base space M is assumed to be a Rie- 
mannian manifold. The isoholonomic problem asks to find the shortest possible piecewise 
smooth loop in M with a given base point Xq G M, that produces a given element Qq of 
the structure group G as its associated holonomy. 

Holonomic structures naturally appear in a mechanical system and have been studied 
from various interests |1]-[H]- Montgomery faced this problem when physical chemists 
attempted to observe the non-Abelian Berry phase (the Wilczek-Zee holonomy) P]-[T^ 
by nuclear magnetic resonance (NMR) experiment. Montgomery fT^ presented various 
formulations of the problem, clarified their relations, and gave partial answers. However, 
even in such an idealized case like a homogeneous bundle, it was difficult to obtain a 
complete solution to the problem, which remained as an open problem to date. 

A decade later after Montgomery's work, the notion of holonomic quantum compu- 
tation was proposed by Zanardi, Rasetti and Pachos [7, "S], in which the Wilczek-Zee 
holonomy is utilized to implement unitary gates necessary to execute a quantum algo- 
rithm. Since then, a large number of researchers ^Sj-llSi have been interested in finding 
control parameters that implement a desired gate. Optimization of the control has been 
an active area of research in view of the decoherence issue. The problem to find the opti- 
mal control is nothing but a typical isoholonomic problem and its solution for an arbitrary 
gate must be urgently provided. 

Let us briefiy review the idea of quantum computation. Quantum computation, 
roughly speaking, consists of three ingredients: (1) an n-qubit register to store infor- 
mation, (2) a unitary matrix U G f/(2") which implements a quantum algorithm, and 
(3) measurements to extract information from the register. In an ordinary implemen- 
tation of a quantum algorithm, we take a system whose Hamiltonian H{X) depends on 
external control parameters A = (A^, . . . , A™). We then properly arrange the parameter 



sequence X(t) as a function of time t so that the desired unitary matrix U is generated as 
a time-evolution operator 



U = Texp 



-'- j\{\{t))dt 



1.11 



where T stands for the time-ordered product. 

Holonomic quantum computing j2], in contrast, makes use of the holonomy associated 
with a loop A(t) in the parameter space. It has been demonstrated ^7] that an arbitrary 
unitary matrix can be implemented as a holonomy by choosing an appropriate loop in 
the parameter space. In fact, there are infinitely many loops that produce a given unitary 
matrix. Here we consider the isoholonomic problem, namely, to find the shortest possible 
loop in the parameter space that yields the given holonomy. This problem has been already 
analyzed previously in Ref . fH] , where various penalty functions useful for numerical search 
for the optimal loop have been employed. Our strategy here is purely geometrical in nature 
and no intense numerical computations are required. In the previous work [T^ we found 
exact optimal loops to produce several unitary gates. In the present paper we extend the 
method of optimal loop construction to implement arbitrary gates. 

This paper is organized as follows. In Section 2 we briefly review the Wilczek-Zee 
holonomy to make this paper self-contained and to establish notation conventions. In 
Section 3 we introduce the geometrical setting for the problem and use it to formulate 
the isoholonomic problem in a variational form. We derive the associated Euler-Lagrange 
equation and solve it explicitly. The solution thus obtained ()3.37|) . which we call the 
horizontal extremal curve, is one of the main results in the first half of the paper. The 
remaining problem is to adjust the solution to satisfy the boundary conditions, namely 
the closed loop condition (|4.H) and the holonomy condition (j4.2j) . This problem is solved 
in Section 4 explicitly and we obtain a set of equations ()4.19|) - ()4.22|) . which we call the 
constructing equations of the controller. These are the main results of this paper and 
are machinery to construct a controller for an arbitrary unitary gate. In Section 5 this 
machinery is applied to several well-known important unitary gates to demonstrate its 
power. Section 6 is devoted to summary and discussions. 

2 Wilczek-Zee holonomy as a unitary gate 

2.1 Wilczek-Zee holonomy 

Here we briefly review the Wilczek-Zee (WZ) holonomy ^U] associated with an adiabatic 
change of the control parameters along a loop in the control manifold. We consider 
a quantum system that has a finite number N of states. Let {H{X)} be a family of 
Hamiltonians parametrized smoothly by A = (A^, . . . , A'") G M, where the set of control 
parameters M is called a control manifold. Eignevalues and eigenstates of H{X) are 



labeled as 

H{X)\l,a;X) = ei{X)\l,a;X) {I = 1, . . . ,L- a = 1, . . . ,ki), (2.1) 

where the l-th eigenvalue £/(A) is fc/-fold degenerate. Assume that no level crossings take 
place, namely, ^^(A) ^ £i'{X) for arbitrary A if / 7^ /'. Then it follows that X]z=i ^i ~ ^■ 
The eigenvectors satisfy the orthonormal condition, {I, a; A|/',/5; A) = 6a'6a/3- It is impor- 
tant to note that there is U{ki) gauge freedom in the choice of {|/, a; A) | a = 1, . . . , ki} at 
each A and /. Namely, we may redefine the eigenvectors by any unitary matrix h G U{ki) 
as 

|Z,«;A)k^5^|/,/3;A)/i^„(A) (2.2) 

/3=1 

without violating the orthonormal condition. 

We adiabatically change the parameters X(t) as a function of time t along a closed 
loop in the control manifold so that A(T) = A(0). It is assumed that the adiabaticity is 
satified, namely, 

{siiXit)) - eiiXit))}T :^ 2nh (2.3) 

is satisfied for / 7^ /' during < t < T. In other words, we change the parameters so 
slowly that no resonant transitions take place between different energy levels |2Uj . 

We will concentrate exclusively on the ground state of the system and drop the index 
/ (= 1) in the following. Accordingly, the basis vectors that span the ground state 
eigenspace are written as |a; A), {a = 1, . . . , k) and arranged in an N x k matrix form as 

\/(A) = (|l;A),|2;A),...,|fc;A)), (2.4) 

which is called an orthonormal /c-frame at A G M. The system evolves, with a given A(t), 
according to the Schrodinger equation 

thj^\Mt)) = H{xmMt))- (2.5) 

Suppose the initial condition is A(0) = Aq and |-?/'a(0)) = \a; Xq). The adiabatic theorem 
PUJ tells us that the state |^a(t)) remains in the ground state eigenspace during the 
time-evolution. Therefore lipait)) is expanded as 

k 

l^<.W) = 5^l/3;A(t))c^,(t). (2.6) 

/3=1 



By substituting ()2.6j) into ()2.5|) . we find 

j^c^ait) = -^£(7(i))c/3a(t) - 5Z ( /5; X{t) 

7=1 



r,X{t))c^a{t), (2.7) 



whose formal solution is 



C/3ait) = exp( — — I e{s)ds ) Texp 



A{s)ds 

/ f3a 



(2.8) 



with the matrix-valued function 



Ap.it) = (d- \{t) ^ a- \{t)^ = X^ ^/?; A 



d 



^\^' 



a\ A 



dX^' 
~~dt' 



(2.9) 



It is easily verified that A*p. 
u(A;)-valued one-form'^ 



-Aap since {\a]\{t))} is orthonormal. We introduce a 






9 



9AA' 



a\\) d\ 



(2.10) 



which is called the Wilczek-Zee (WZ) connection. Then the unitary matrix appearing in 
(12. 8 j) is rewritten as 



r{t) =Vexp 



m 



A], 



f2.ir 



A(0) 



where V stands for the path-ordered product. As noted in ()2.2|) the frame ()2.4|) can be 
redefined by a family of unitary matrices h{X) G U{k). The WZ connection transforms 
under the change of frame as 



(2.12) 



A^A' = h^Ah + h^dh. 



This is nothing but the gauge transformation rule of a non-Abelian gauge potential 

We assumed that the control parameter \{t) comes back to the initial point A(T) = 
A(0) = Aq. However, the state \'il)a{T)) fails to assume the initial state and is subject to 
a unitary rotation as 



IV'.(T)) =exp(-- / e{s)ds]Y,\^M)rp.{T). 



The unitary matrix 



r[\] := r{T) = Pexp (-([ a) e U{k) 



(2.13) 



(2.14) 



is called the holonomy matrix associated with the loop A(t). It is important to realize 
that r[X\ is independent of the parametrization of the loop A(t), namely, it is independent 
of how fast the loop A is traversed, so long as the adiabaticity is observed, and that it 
depends only on the geometrical image of A in M. 

^We denote the Lie algebra of the Lie group U{k) by u{k), which is the set of fc-dimensional skew- 
Hermite matrices. 



2.2 Quantum computation with holonomy 

In quantum computation one implements a quantum algorithm by a product of various 
unitary gates. It is a natural idea to use the WZ holonomy to produce unitary gates 
necessary for quantum computation. Zanardi and Rasetti [2] were the first who proposed 
this holonomic quantum computation (HQC). To implement an n-qubit resistor we take 
a quantum system whose ground state is fc-fold degenerate where k = 2". We call the 
A^- dimensional Hilbert space a working space and call the /c-dimensional subspace a qubit 
space. Then by changing the control parameter adiabatically we will obtain any unitary 
gate as a resultant holonomy ()2.14|) . Of course we need to design an appropriate control 
loop A to implement a particular unitary gate. It is easy, in principle, to compute the 
holonomy for a given loop. In contrast, to find a loop A which produces a specified unitary 
matrix F as its holonomy is far from trivial. Moreover, to build a working quantum 
computer it is strongly desired to reduce the time required to manipulate the computer 
since a sequence of operations should be carried out before decoherence extinguishes 
quantum information from the system. At the same time, the control parameter must 
be changed as slowly as possible to keep adiabaticity intact. Therefore our task is to 
find a control loop as short as possible to fulfill these seemingly opposed conditions. This 
is a typical example of the so-called isoholonomic problem, which is first formulated by 
Montgomery jj. In the next section we introduce a geometric setting in a form suitable 
for our expositions. 

3 Formulation of the problem and its solution 

3.1 Geometrical setting 



The WZ connection is identified with the canonical connection |22j of the homogenous 
bundle, as pointed out by Fujii [221 • While precise definitions of these terms can be found 
in books fn\ 1^ . we outline the geometrical setting of the problem here to make this 
paper self-contained. 

Suppose that the system has a family of Hamiltonians acting on the Hilbert space C^ 
and that the ground state of each Hamiltonian is k-iold degenerate {k < N). The most 
natural mathematical setting to describe this system is the principal bundle (S^^kiC), 
GN,k{^),T^,U{k)), which consists of the Stiefel manifold Sjy^ki^), the Grassmann mani- 
folds GN,k{'C), the projection map n : SN,k{'C) -^ GN,k{C), and the unitary group U{k) 
as explained below. 

The Stiefel manifold is the set of orthonormal fc-frames in C^, 

SnAQ = {V^ M{N, k- C) I VW = 4}, (3.1) 

where M{N, k; C) is the set oi N x k complex matrices and I^ is the /c-dimensional unit 



matrix. The unitary group U{k) acts on SN,k{'C) from the right 

SmAQ X U{k) ^ S^,,(C), {V, h) ^ Vh (3.2) 

by means of matrix product. It should be noted that this action is free. In other words, 
h = Ik ii there exists a point V G SN,k{C) such that Vh = V. 

The Grassmann manifold is defined as the set of fc-dimensional hyperplanes in C^, 

GnAQ = {P^ M{N, N; C) I P^ = P, P^ = P, trP = k}, (3.3) 

where P is a projection operator to a hyperplane in C^ and the condition trP = k 
guarantees that the hyperplane is indeed /c-dimensional. 
The projection map vr : SN,k{'C) — ^ GN,k{C) is defined as 

TT : 1/ ^ P := VVl (3.4) 

It is easily proved that the map tt is surjective. Namely, for any P G GN,k{C), there is 
V G SN,k{C) such that niV) = P. The right action of /i G U{k) sends a point V G S'Ar,fc(C) 
to a point V/i on the same fiber since 

7r{Vh) = {Vh)(yhy = Vhh^V^ = VV^ = TciV). (3.5) 

Thus the Stiefel manifold SN,k{'C) becomes a principal bundle over GN,k{C) with the 
structure group U{k). 

Moreover, the group U{N) acts on both SN,k{C) and GN,k{C) as 



U{N) X SnAQ - 


■* S'Af,fc(C), 


(^, ^) ^ ^v^, 


(3.6) 


UiN) X Gjv,fe(C) - 


-^ G'Af,A:(C), 


((7,P)H^(7P(7t 


(3.7) 



by matrix product. It is easily verified that TT{gV) = g7r{V)g^. This action is transitive, 
namely, there is (? G U{N) for any V,V' G SnA^) ^^ch that V = gV. There is also 
g G U{N) for any P, P' G GN,k{C) such that P' = gPg^. The stabilizer group of each point 
in SN,k{'C) is isomorphic to U{N — k) while that of each point in GnA'^) i^ isomorphic 
to U{k) X U{N — k). Thus, they are homogeneous spaces and the fiber bundle 

TT : SnAQ = U{N)/U{N -k)^ GnA^) = U{N)/{U{k) x U{N - k)) (3.8) 

is call a homogeneous bundle. 

The canonical connection form on S^A'^) i^ defined as a u(/i;)-valued one- form 

A = V^dV, (3.9) 

which is a generalization of the WZ connection ()2.10|) . This is characterized as the unique 
connection that is invariant under the action ()3.6p . The associated curvature two- form is 
then defined as 

F = dA + AAA = dV^AdV + VUV A VUV = dV^ A {In - VV^)dV. (3.10) 



These manifolds are equipped with Rieinannian metrics. We define a metric 

\\dVf = tT{dVUV) (3.11) 

for the Stiefel manifold and 

\\dPf = ti (dPdP) (3.12) 

for the Grassmann manifold. 

3.2 The isoholonomic problem 

Here we reformulate the WZ holonomy in terms of the geometric terminology introduced 
above. The state vector ip(t) G C^ evolves according to the Schrodinger equation 

^hJ^m = H{t)m- (3.13) 

The Hamiltonian admits a spectral decomposition 

L 

H{t) = Y,ei{t)Pi{t) (3.14) 

1=1 

with projection operators Pi(t). Therefore, the set of energy eigenvalues (ei, . . . ,6^) and 
orthogonal projectors (Pi,...,Pj.) constitutes a complete set of control parameters of 
the system. Now we concentrate on the eigenspace associated with the lowest energy, 
which is assumed to be identically zero, ei = 0. We write -Pi(t) as P{t) for simplicity. 
Suppose that the degree of degeneracy k = ti P{t) is constant. For each t, there exists 
V{t) G SN,k{C) such that P(t) = V{t)V'^{t). By adiabatic approximation we mean 
substitution of ip{t) G C^ by a reduced state vector (f)(t) G C'^ as 

m = V{t)<P{t). (3.15) 

Since H{t)ip{t) = eiip{t) = 0, the Schrodinger equation ()3.13|) becomes 

f...>)^o (3.ia) 

and its formal solution is written as 

(l){t) = Vexp(- fvUVjcpiO). (3.17) 

Therefore ^(t) is written as 

^(t) = V{t)Vexp(- fvUVj V\0)ij{0). (3.18) 



In particular, when the control parameter conies back to the initial point as P{T) = P(0), 
the holonomy F G U{k) is defined via*^ 

^{T) = V{Q)rm (3.19) 

and it is given explicitly as 

r = V{OyV{T)Vexp(- fvUVj . (3.20) 

If the condition 

V^^ = (3.21) 

dt ^ ' 

is satisfied, the curve V{t) in SN^ki^C) is called a horizontal lift of the curve P{t) = n{V{t)) 

in GN,k{'C). Then the holonomy ()3.2()|1 is reduced to 

r = V^{0)V{T) G U{k). (3.22) 

Now we are ready to state the isoholonomic problem in the present context; given a 
specified unitary gate ?7gate £ U{k) and a fixed point Pq G GN,k{'C), find the shortest 
loop P{t) in Gjv^fc(C) with the base points -P(O) = P{T) = Pq whose horizontal lift V{t) 
in Sisi^kiC) produces a holonomy P that coincides with t/gate- This problem was first 
motivated from experimental study of geometric phase and investigated in detail from a 
mathematician's viewpoint by Montgomery pp. 

We now formulate the isoholonomic problem as a variational problem. The length of 
the horizontal curve V{t) is evaluated by the functional 

where i7(t) G u{k) is a Lagrange multiplier to impose the horizontal condition 1)3.211) on 
the curve V{t). Note that the value of the functional S is equal to the length of the 
projected curve P(t) = Tr(y(t)), 

S=riJ'-f'-^),L (3.24) 

Jo 2 \dt dt J ^ ' 

Thus the problem is formulated as follows; find a curve V{t) that attains an extremal 
value of the functional ()3.23|) and satisfies the boundary condition ()3.22j) . 

'^The definition of the liolonomy presented here is shghtly different from the one given in the previous 
Letter (19, . To make a correct sense as a miitary gate the holonomy is to be defined in the present form. 



3.3 The solution: horizontal extremal curve 

Our task is to find a solution of tlie variational problem of the functional ()3.23p . Now 
we derive the associated Euler-Lagrange equation and solve it explicitly. A variation 
of the curve V{t) is defined by an arbitrary smooth function r]{t) G u{N) such that 
ri{0) = T]{T) = and an infinitesimal parameter e G M as 

V,{t) = {l + ev{t))Vit). (3.25) 

By substituting K(t) into (J3.23J) and differentiating with respect to e, the extremal con- 
dition yields 

e=0 



0=^ 

de 



r -,t=T 

= tr {r](VV^ - VV^ - VQV^)} 
L Ji=o 

tr {vjiVV^ - VV^ - VnV^) } dt. (3.26) 

Thus we obtain the Euler-Lagrange equation 

^{VV'^ -VV'^ \VnV^) = ^. (3.27) 

We reproduce the horizontal equation V'^V = from the extremal condition with re- 
spect to i7(t). Finally, the isoholonomic problem is reduced to the set of equations 
()3.21|) and ()3.27|) . which we call a horizontal extremal equation. It may be regarded 
as a homogeneous-space version of the Wong equation |24j . 

Next, we solve the equations (j3.21|) and (J3.27J) . The equation (j3.27p is integrated to 
yield 

VV^ - VV^ + VQV^ = const = X G u(A^). (3.28) 

Conjugation of the horizontal condition ()3.21|) yields V^V = 0. Then, by multiplying V 
on ()3.28|) from the right we obtain 

V + VQ = XV. (3.29) 

By multiplying V'^ on (J3.29J) from the left we obtain 

Q = V^XV. (3.30) 

The equation ()3.29p implies V = XV — Vf2, and hence the time derivative of f2{t) becomes 

Q = V^XV + V^XV = V^X{XV - VQ) + (-1/^X + f2V^)XV = [Q, Q] = 0. (3.31) 

Therefore, f](t) is actually a constant. Thus the solution of ()3.29|) and ()3.30|) is 

V{t) = e*^ Vo e"*^, Q = V^XVo. (3.32) 
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We call this solution the horizontal extremal curve. Then ()3.28|) becomes 

{XV - vf2)v^ - v{-v^x + ov^) + vnv^ = x, 

which is arranged as 

X -{VV^X + XVV^ -VV^XVV^)=0, (3.33) 

where we used fl3.30|) . We may take, without loss of generality, 

Vo=( ^']e SnAQ (3.34) 



^ 
as the initial point. We can parametrize X G u(A^), which satisfies fl3.30|) . as 

X ^ ( _« , ^; ) (3.35) 

with W G M{k, N — k]C) and Z G u(A^ — k). Then the constraint equation ()3.33|) forces 
us to choose 

Z = 0. (3.36) 

Finally, we obtained a complete set of solution (J3.32J) of the horizontal extremal equa- 
tion (j3.21|) and (J3.27J) . When we take the initial point Vq as (J3.34J) . the solutions are 
parametrized by constant matrices i7 G u{k) and W G M{k,N — k;C). For definiteness 
we write down the complete solution 

V{t) = e^^Voe-^^, X=^_^, ^). (3.37) 

This is one of our main results. We call the matrix X a controller. At this time the 
holonomy ()3.22|) is expressed as 

r = V^{0)V{T) = V^ e^^ Vo e-™ G U{k). (3.38) 

These results ()3.37p and ()3.38|) have been also given in Montgomery's paper*^. In the 
present paper we took a different approach from his. Here we wrote down the Euler- 
Lagrange equation and solved it directly. 

We evaluate the length of the extremal curve for later convenience by substituting 
(HOTIl into (IT^ as 



'^In his paper ^ Montgomery cited Bar's theorem to complete the proof. However, Bar's paper being 
a diploma thesis, it is not widely available. Therefore we took a more direct approach to justify them. 
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4 Solution to the inverse problem 

Once the solution (|3.37|) of the horizontal extremal equation is obtained, the remaining 
problem is to find the matrices i7 and W that satisfy the closed loop condition 

ViT)V\T) = e^^V^oKi^e-^^ = VoV^ (4.1) 

and the holonomy condition 

V^ y (T) = V^ e^^ Vo e-^^ = [/gate (4.2) 

for a specific unitary gate t/gate ^ U{k). Montgomery P presented this inverse problem 
as an open problem. In this section we give a scheme to construct systematically a series 
of solutions to this problem and in the next section we will apply it to implement various 
important unitary gates. 

4.1 Equivalence class 

There is a class of equivalent solutions with a given initial condition Vq and a given 
final condition V{T) = Vof/gatc- Here we clarify the equivalence relation among solutions 
{V{t)} that have the form (HnTfll and satisfy dHH) and (Q. 

We say that two solutions V{t) and V'{t) are equivalent if there are elements g G U{N) 
and h G U{k) such that V{t) and 

V'it) = gVit)hl (4.3) 

satisfy the same boundary conditions 

gVoh^ = Vo (4.4) 

and 

hUg^teh^ = t/gatc- (4.5) 

For the initial point ()3.34|) . the condition ()4.4j) states that g G U{N) must have a block- 
diagonal form 

■ hi 



^ , ,, h = hieU{k), h2eU{N-k). (4.6) 

U ri2 
The controller X' of V'{t) are then found from 

V'{t) = gV{t)h^ = ge'^g^gVoh^he~'^h^ = e'^^^^ gVoh^e-'^^''^ = e'^^^Vo e-*^^'^\ (4.7) 

In summary, two controllers X and X' are equivalent if and only if there are unitary 
matrices hi G U{k) and /i2 G U{N — k) such that 



x-( ^ iy\ ,_( hiQhl hiWh\ \ ,jj,^_jj 



(4.8) 
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4.2 ^(1) holonomy 

Here we calculate the holonomy for the case N = 2 and k = 1. In this case the homoge- 
neous bundle vr : 5*2, i(C) — > G2,i(C) is the Hopf bundle n : S^ -^ S'^ with the structure 
group U{1) and the WZ holonomy reduces to the Berry phase. In the subsequent sub- 
section we will generalize this result to a non-Abelian holonomy. We normalize the cycle 
time as T = 1 in the following. Using real numbers ^1,^2,^3 G M we parametrize the 
controller as 

^, / 2iW3 iWl + W2 \ . r ■ ■ ■ / A r^\ 

X = \ . ^ = IW^l + IWiCJi + lW2(y2 + «W30"3, (4.9) 

\^ IWi — W2 J 

where {cj} are the Pauli matrices. Its exponentiation is 

gix ^ git«,3 (^ J cos pt + in-cr sin pt) , (4. 10) 

where p and n are defined as 

p := ||iu|| = ^y{wly + (^2)^ + (ws)^, it> = ||if ||n. (4.11) 

The associated horizontal extremal curve ()3.32j) then becomes 

V{t) = e*^ Vo e-^ = e--3 ( ^ospt + ^n,smpt \ 
^ ' \ [tn I — n2) sm pt J ^ ' 

and the projected curve in S'^ becomes 

Pit) = vm^t) 

= -I + -cr ■ [n{n ■ 63) + (63 — n{n ■ 63)) cos 2pt — (n x 63) sin 2pt] , (4.13) 

where 63 = (0, 0, 1). We see from ()4.13p that the point P{t) in S'^ starts at the north pole 
63 of the sphere and moves along a small circle with the axis n in the clockwise sense by 
the angle 2pt. The point P{t) comes back to the north pole when t satisfies 2pt = 2iTn 
with an integer n. To make a closed loop, namely, to satisfy the loop condition ()4.1|) at 
t = T = 1, the control parameters must satisfy 

p=\\w\\=mv (n = ±1,±2,...). (4.14) 

Then, the point P(t) travels the same small circle n times during < t < 1. Therefore, 
the integer n counts the winding number of the loop. At t = 1, cosp = (—1)" and the 
holonomy ()4.2|) is evaluated as 

V^o^ e^ Vo e-^ = e'""^ {-!)'' e-^'""^ = e'^^^^""'^) = f/g^te = e'^. (4.15) 

Thus, to generate the holonomy Legate = e*''', the controller parameters are fixed as 



W3 = nvr — 7, wi + iw2 = e ^'^\/{mi)'^ — (mr — 7)2. (4.16) 
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This is the solution to the inverse problem defined by ()4.H) and ()4.2|) . Here the nonvanish- 
ing integer n must satisfy (ravr)^ — (nvr — 7)^ > 0. The real parameter is not fixed by the 
loop condition and the holonomy condition. The phase /12 = e**^ parametrizes solutions 
in an equivalence class as observed in ()4.8|) . The integer n classifies inequivalent classes. 
The length of the loop, (|3.39|) . is now evaluated as 

S = ii {W^W)T = {rn^f -{n-K- -if . (4.17) 

For a fixed 7 in the range < 7 < 27r, the simple loop with n = 1 is the shortest one 
among the extremal loops. Thus, we conclude that the controller of f/gate = e*'^ is 



V ie-'^^T:^ - (tt - 7)2 y ^ ' 

We call this solution a small circle solution because of its geometric picture mentioned 
above. 

4.3 U{k) holonomy 

Here we give a prescription to construct a controller matrix X that generates a specific 
unitary gate Legato- In other words, we give a systematic method to solve the inverse 
problem (j4.2j) . It turns out that the working space should have a dimension N > 2k 
to apply our method. In the following we assume that N = 2k. The time interval is 
normalized as T = 1 as before. 

Our method consists of three steps: first, diagonalize the unitary matrix Legate to be 
implemented, second, construct a diagonal controller matrix by combining small circle 
solutions, third, undo diagonalization of the controller. 

In the first step, we diagonalize a given unitary matrix Ug^tc £ U{k) as 

i?^f/gatei? = t/diag = diag(e^^S . . . , 6*'^'=) (4.19) 

with R G U{k). Each eigenvalue 7^ is taken in the range < 7^ < 2-7?. In the second step, 
we combine single loop solutions associated with the Berry phase to construct two k x k 
matrices 

^diag = diag(ia;i, . . . , iuk), ujj = 2{tt - 7^), (4.20) 



Wdiag = diag(2ri, . . . , ^Tfc), r, = e'^^ ^tt^ - (tt - ^,'f. (4.21) 

Then we obtain a diagonal controller 



-^diag - I _yT.t 



diag ''^diag 



^diag 
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In the third step, we construct the controller X as 

which is a 2 A; X 2k matrix. We call the set of equations, (pTT9|l . (jODI), ^4~T\} and (gSS), 
constructing equations of the controller. This is the main result of this paper. 

It is easily verified that the controller X constructed above satisfies the holonomy 
condition ()4.2|) . The diagonal controller Xjiag is actually a direct sum of controllers 
()4.18|) . which generate Berry phases {e*'''^}. Hence, its holonomy is also a direct sum of 
the Berry phases (j4.15p as 

and hence we have 



V^ e^ Vo e~^ = RV^ e^^-^ Vq R^Re-^^'^^R^ = RU^i^^R^ = U^ 



■ate- 



5 Optimizing holonomic quantum computation 

Now we apply the prescription developed so far to construct controllers of several specific 
unitary gates, which are fundamental ingredients of quantum computation. Our examples 
are the Hadamard gate, the CNOT gate, and the two-qubit discrete Fourier transformation 
(DFT) gate. For each unitary gate t/gate, we need to calculate the diagonalizing matrix 
R. Then the constructing equations of the controller, ()4.19|) - ()4.22p . provide the desired 
optimal controller matrices. 

5.1 Hadamard gate 

The Hadamard gate is a one-qubit gate defined as 



It is diagonalized by 



cos ^ ^ sm g 



sm ^ cos g 



R= .._^ ___ (5-2) 



as 



Needless to say. 



R^Uu^^R = r J _M . (5.3) 

71 V2 + V2 . n V2-V2 

cos— = , sm — = . (5.4) 

8 2 8 2 ^ ^ 

Therefore, we have 71 = and 72 = vr. We may put 0i = 02 = 0. The ingredients of the 
constructing equations of the controller, ()4.19|) - ()4.22|) . are calculated as 

^diag = diag(227r, 0), IVdiag = diag(0,i7r), (5.5) 
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and hence 



RildiapR — — = 



v^+1 1 



tTC 



RW.^.„ = -r- 



-V2-v^ 



Substituting these into ()4.22|) . we obtain the optimal controller of the Hadamard gate. 
5.2 CNOT gate 

One of the most important 2-qubit gates is the CNOT gate defined as 



(5.6) 



U. 



CNOT 



It is diagonalized by 



V2 



/ 1 \ 

10 

1 

\ 1 / 



/ ^2 \ 

v^ 

1-1 

V 11/ 



(5.7) 



(5.^ 



as 



R^U, 



CNOT 



R 



/ 1 \ 

10 

10 

\ -1 / 

Therefore, we have 7i = 72 = 73 = and 74 = vr. The ingredients of the controller are 



(5.9) 



i^diag = diag(2m, 2i'K, 2iTT, 0), VFdiag = diag(0, 0, 0, itt), 

/ \ 



(5.10) 



Rf2d\f,tyR 



tTC 



(5.11) 



and hence 

/ 2 \ 

2 _ iir 

11 ' ^'"s ~ 72 0-1 

\0011/ \0001/ 

Substituting these into ()4.22|) . we obtain the optimal controller of the CNOT gate. 

5.3 DFT2 gate 

Discrete Fourier transformation (DFT) gates are important in many quantum algorithms 
including Shor's algorithm for integer factorization. The two-qubit DFT (DFT2) is a 
unitary transformation 



U^ 



I 



DFT2 



/I 1 1 1 \ 

1 i -1 -i 

1-11-1 

\ 1 -i -1 i / 



(5.12) 
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It is diagonalized by 



R 



/ 1 v^ -1 \ 

1 1 -V2 

-1 V2 1 

V 1 1 V2 J 



(5.13) 



as 



R ?7dft2-R 



(5.14) 



/ 1 \ 

10 

0-10 

\ i J 

Therefore, we have 7i = 72 = 0, 73 = vr and 74 = tt/2. Thus the ingredients of the 
controller are 



^diag = diag(2i7r, 2i7T, 0, ivr), PVdiag = diag(0, 0, ivr, i7rV3/2), 



(5.15) 



and hence 



RQdmsrR — — 



^diag-' 



/ 3 1 1 1 \ 
12-10 
1-13-1 

\ 1 -1 2 / 



777- 



/ -1 \ 

1 -v/372 

1 

Vo 1 v^ I 



(5.16) 
Substituting these into ()4.22j) . we finally obtain the optimal controller of the DFT2 gate. 



6 Summary and discussions 

Let us summarize our argument. We briefiy reviewed the WZ holonomy and discussed that 
it may be utilizable for implementation of quantum computation. The WZ holonomy is 
neatly described in terms of differential geometry of a homogeneous bundle, which consists 
of Stiefel and Grassmann manifolds and is equipped with the canonical connection. We 
formulated the optimization problem of control in holonomic quantum computation in a 
form of the isoholonomic problem in the homogenous bundle. We would like to emphasize 
that it had been left unsolved for more than a decade after the first proposal. We derived 
a set of equations, ()3.2H) and ()3.27p . that characterizes the optimal control and solved 
it to obtain the horizontal extremal curve ()3.37|) . The curve must satisfy two boundary 
conditions, (J4.H) and (j4.2j) . to be a closed loop in the control manifold and to produce a 
specified unitary gate as a holonomy. We solved this inverse problem by combining small 
circle solutions ()4.18p to U{1) holonomy into a direct sum. We provided a prescription 
fl4.19|) - ()4.22|) to construct exactly an optimal controller for any unitary gate. Finally we 
applied our prescription to several important quantum gates, transform gate. 

We would like to discuss prospective development of the results presented above. Al- 
though our prescription is applicable to arbitrarily large qubit gates, the homogeneous 



17 



bundle seems rather over-idealized for practical applications. A realistic quantum system 
may have smaller control manifold M than the Grassmann manifold. The restricted 
control manifold M is embedded in the Grassmann manifold by an embedding map 
f : M —^ GN,k{'C) and we need to study the isoholonomic problem in the puUbacked 
bundle /*S'7v,fc(C). Furthermore, the available working Hilbert space in a realistic system 
may not have dimensions as large as A^ > 2k. Actually, even when A^ < 2k, sequential 
operations of single loop solutions can generate any unitary gate. However, such a patched 
solution could not be optimal. These problems will be treated separately in our future 
publications. 
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